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ABSTRACT 

The current models of early star and galaxy formation are based upon the 
hierarchical growth of dark matter halos, within which the baryons condense 
into stars after cooling down from a hot diffuse phase. The latter is replenished 
by infall of outer gas into the halo potential wells; this includes a fraction 
previously expelled and preheated, due to momentum and energy fed back 
by the SNe which follow the star formation. We identify such an implied hot 
phase with the medium known to radiate powerful X-rays in clusters and in 
groups of galaxies. We show that the amount of the hot component required by 
the current star formation models is enough to be observable out to redshifts 
z ~ 1.5 in forthcoming deep surveys from Chandra and XMM, especially in case 
the star formation rate is high at such and earlier z. These X-ray emissions 
constitute a necessary counterpart, and will provide a much wanted probe of the 
SF process itself (in particular, of the SN feedback), to parallel and complement 
the currently debated data from optical and IR observations of the young stars. 



Subject headings: galaxies: formation - galaxies: clusters: general - intergalactic 
medium - X-rays: general 
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1. Introduction 

Is there any link of the early star formation rate (SFR) with extended extragalactic 
sources expected to appear in the deep X-ray surveys planned (see Giacconi 1998) with 
Chandra and with XMM 

We expect some such connection to exist, if we carry a step further the views started 
by several teams (Munich: Kauffmann, White & Guiderdoni 1993; Durham: Cole et al. 
1994, Baugh et al. 1998; Santa Cruz: Somerville & Primack 1998) to describe the processes 
that lead to galaxy and star formation. These are based upon the hierarchical growth of 
structures gravitationally dominated by dark matter (DM) halos; from the size of a galaxy 
to that of a galaxy group and then of a cluster, the growth occurs by repeated merging of 
smaller into larger structures. These views envision the baryons as condensing into stars at 
the minima of the forming DM potential wells, upon cooling down from a diffuse hot phase 
at the virial temperature of the wells. 

We will show that such implied hot component - constituting, in fact, the intracluster 
medium (ICM) in the larger halos - yields copious bremsstrahlung emissions observable in 
X-rays. These are closely linked with the optical or infrared stellar light from which the 
SFR is currently inferred. The link is provided by the energy and momentum fed back into 
the hot phase by the stellar winds and the supernovae (SNe) following the star formation. 
This process is especially relevant in the early, small but dense halos; some feedback is 
essential there to prevent the runaway cooling of all the contained baryons, the so-called 
cooling catastrophe (White & Rees 1978). 

2. Key quantities 

SN explosions, with some contribution from stellar winds (Bressan, Chiosi & Fagotto 
1994), provide an energy output E* = E S n Vsn Am*; here E SN 10 51 erg is the energy of 
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a Type II SN explosion, and r] SN fx 4 10 3 Am* is the combined efficiency for making SNe 
and massive blue stars when the mass Am* condenses into the IMF. 

Two key parameters will gauge the effectiveness of the feedback against the depth of 
the containing potential well. Recall that the depth, for a halo of mass M with internal 
density p oc (1 + z) 3 following the background's, is marked by the circular velocity 
v = yj GMj R v oc M 1 / 3 p 1 / 6 , or by the virial temperature kT v oc M/R v oc v 2 , which takes on 
values around 4 keV in rich clusters. 

One parameter will gauge the dynamical effect of stars and of SN explosions by way of 
momentum transfer onto the surrounding gas; the resulting galactic winds can directly expel 
a fraction Ara^/ra^ = /*, especially from shallow potential wells. The proper parameter is 
the fractional energy converted into bulk kinetic energy at the escape velocity 

e = f*m h v 2 /E* . (1) 

This will take on values of order 10 _1 in wells with v ~ 150 km/s; but equally important 
will be its differential behavior in shallower and deeper wells. 

The other parameter T v /T* will gauge the importance of the thermal effect of the SNe. 
Here kT v measures the gravitational energy of the outer baryons, amenable to conversion 
into heat as they infall into the wells. The "stellar" temperature kT* = (1 — e )E ie m p /3 Am h 
measures the preheating level of the expelled fraction, provided by the SNe. This may be 
preliminarly estimated as 

kT, = 0.7 (1 - e Q ) Am*/Am h £ 0.2 keV , (2) 

for a stellar baryonic fraction exceeding 1/5, the value appropriate only for rich clusters (see 
Renzini 1997). But kT v fx 0.2 keV is also the virial value corresponding to M ~ 5 10 12 M ; 
so we may surmise that somewhere in the range from a poor group to a large galaxy stellar 
preheating starts to prevent the gas from piling up within the wells. Such masses begin to 
virialize at z ~ 2 ± 0.5 in current hierarchical cosmogonies. 
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We shall investigate whether more SF at early z, which requires more baryons 
condensed within the small halos, also leads to more hot gas retained and so to stronger 
and widespread X-ray emissions. To proceed beyond estimates we need a model for the 
partition, shifting with M and z, of the baryons among the condensed, the cool, the expelled 
and the hot components. 

3. The model 
3.1. Optical emissions from stars 

To that effect we start from the semi-analytic models (SAMs) developed by the teams 
of Munich, Durham and Santa Cruz. These include the "merging histories" of DM halos 
as they grow hierarchically through stochastic merging events; the above authors describe 
such histories using Monte Carlo simulations. We use instead the equivalent analytic 
probabilities provided by Bower (1991) and by Lacey & Cole (1993), the so-called extended 
Press & Schechter theory. This is because our model is considerably more complex as it 
covers also the X-ray emitting baryons; so we need to cut down the computer time required 
to average our many observables by convolutions over the merging histories. 

Along with such a rendition of the DM dynamics, the numerical package devised and 
built by one of us (N. Menci) includes the current SAM "recipes" to follow the condensed, 
the cool and the hot baryons. Here we adopt the Durham implementation where: the mass 
of cool gas Am c is what can cool down at the center of the halos over their survival time; 
the SF from the cooling stuff obeys Am* = [m c — m* — m h ] At/r° (v/300 bis" 1 )*; the 
fraction returned from the cool to the hot component by galactic winds is given by Am), 
= Am*(vh/v) ah , corresponding to an effective e c = Am h v 2 /E* oc v 2 ~ ah . The luminosities of 
the stars so produced are computed on convolving m* with the spectral energy distributions 
provided by the model of stellar population synthesis by Bruzual & Chariot (1998). 

The primary parameters appearing here are t *, a*, a^, Vh', reasonable values for these 
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are known to produce reasonable first approximations to the galaxy Tully-Fisher relation, 
and to their colors and related luminosity functions. In closer detail, the latter are improved 
(as discussed in Cole et al. 1994; Somerville & Primack 1998) on introducing secondary 
parameters to modulate the shape of the bright end (like the amount of absorbing dust) 
and the shape of the local, faint end (like the coalescence rate of the galactic baryonic 
cores within the halos). Such galactic parameters are implemented in our package after the 
Durham version, but in fact do not affect the hot diffuse ICM pervading the halos. 

In this Letter we spare details and figures concerning the O and IR observables used 
here as a calibration, and refer to Menci & Cavaliere (1999; MC99) and Poli et al. (1999). 
We confine ourselves to show in fig. 1 (top) our results concerning the intrinsic SFR 
corresponding to prompt blue stellar light. When we use the original Durham values for the 
parameters, we obtain optical results equal to theirs, see solid line in fig. 1 (top) and the 
related caption. 



3.2. X-ray emissions from the hot phase 

Our new step is the X-ray emission from the hot component (ICM) required by such 
SFRs. In terms of the ICM temperature T and particle density n, the continuum luminosity 
reads L x oc n 2 R x T 1 / 2 . Using T T v oc p-^M 1 / 3 the scaling is conveniently recast as 

L x oc p 1/2 (m h /M) 2 (n 2 /m) 2 J T 2 . (3) 

The shape factor J depends only weakly on z and on halo mass M; it arises from the 
volume integration over the inner ICM profiles n 2 (r) and T^ 2 (r), computed (see CMT99) 
on assuming for the ICM hydrostatic equilibrium in the DM potential provided by Navarro, 
Frenk & White (1997). The main dependences arise from expulsion and heating that 
affect the fraction mh/M of hot ICM inside the wells, and its internal density n 2 relative 
to the external n\. In fact, taking constant n 2 /ni and rrih/M would give Lx oc T 2 , the 
scale-invariant relation for ICM just following the gravitational, hierarchical clustering; 
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but the result would clash against the observations referred to below. To break such 
scale-invariance the ICM must respond actively under the drive of the stellar feedback. 

The thermal feedback affects the ratio n 2 /ni, which depends on the strength of the 
accretion shocks (see CMT97) induced at about the virial radius R v when the outer gas 
falls supersonically into the halo M. The latter is built up hierachically by merging with, 
or by accretion of lumps of mass M'; the associated gas is at a temperature T' . Given this, 
the density ratio across the shock depends on its strength T/T' after the functional form 
g(T/T') given by CMT99; this rises steeply from 1 when T/T' is moderately increased, but 
saturates to 4 for large values of T/T'. 

The dynamical feedback modulates this behavior, as it yields two values for T'\ the 
"stellar" value T[ applies to the fraction /* = Am h /mh ejected beyond the virial radius of 
M' and preheated; the gravitational applies to the complementary 1 — /* that remains 
virialized inside the lump. Both such components fall into the main well, and the weighted 
density ratio is given by (MC99) 

K/ni) 2 = /, g\T/T[) + (1 - /*) g\T/T' v ) . (4) 

Meanwhile, some adiabatic compression takes place in the settling of the shocked ICM 
into the well. Compressions and shocks occur repeatedly during the hierarchical growth. 
Their balance may be expressed with a polytropic equilibrium, implying the radial profile 
T(r) oc n r ' 1 (r) with 1 < T < 1.3 (see Cavaliere & Fusco-Femiano 1978, Tozzi & Norman 
1999); this converts the previous history of shocks and compressions in the cluster into a 
space stratification. Such ICM profiles are used to compute the shape factor T. 

Averaging eq. (4) over all merging histories (i.e., over the distributions of M') that 
lead to M yields the bent L x — T relations in fig. 1 (middle); these start steep at the scale 
of groups, gradually flatten, and saturate toward L x oc T 2 only in very rich clusters, in 
agreement with the observations by Ponman et al. (1996); Mushotzky & Scharf (1997); 
Markevitch (1998); Allen & Fabian (1998). Such a bent shape arises basically because in 
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groups the stellar preheating temperature T[ matches or exceeds the virial T as anticipated 
in Sect. 2; then T/T' and n 2 /nx approach 1, so L x oc n^/nf is weak; in fact, it is weaker 
when the expelled/heated fraction is larger. In rich clusters instead, T exceeds any 
external and the shocks are always strong; there n^jrix saturates to 4, and Lx saturates 
toward L x oc T 2 . Additionally, a larger expelled fraction /* decreases the relative amount 
m h /M of ICM inside the shallow wells, and again weakens L x in groups. 

4. Results 

All that joins nicely with the SAMs, which are based on the same hierarchical merging 
histories as said above. In fact, our treatment of the hot ICM can be grafted onto the 
SAMs because it proceeds from the same basic equations and uses the same parameters. 
A double bonus of our approach is that heating the gas while expelling it outwards to low 
densities is a process best suited to generate in groups the high "entropy floor" observed by 
Ponman Cannon & Navarro (1999); in fact, we obtain kT,/n 2 / 3 ~ 0.5 10 8/3 ~ 10 2 keV cm 2 . 
Such a value is conserved near the center of forming clusters by the subsequent adiabatic 
compressions. But as a cluster forms, increasingly strong shocks develop farther out and 
deposit the entropy AS* oc InT n 2 / 3 /T'n^ 3 ; so the entropy rises outwards, in accord with 
what is observed in clusters (David, Jones & Forman 1996; Ponman et al. 1999). 

The derived Lx — T relation constitutes the intermediate step (analogous to the 
Tully-Fisher relation for the star light) to compute the X-ray luminosity functions and the 
source counts using standard formulae. We have included the line emissions using a routine 
updating Raymond & Smith (1977), kindly provided by P. Tozzi. 

Our fiducial cosmogony/cosmology is provided by standard cold DM initial 
perturbations in a flat, low-density Universe with Q Q = 0.3 and Q\ = 0.7, and H =70 
km/s Mpc. Not only this is indicated by the data concerning SNe la, but also it makes the 
baryonic fraction il b = 0.04 suggested by the cosmological nucleosynthesis (see Olive 1998) 
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consistent with the values flb/^o around 0.15 from mass inventories in clusters (White et 
al. 1993; Ettori & Fabian 1999). 

So equipped, we explore the counterparts in X-rays of two extreme evolutions of the 
SFR for z > 1 (shown in fig. 1, top), which bracket the currently debated shapes derived 
from the O-IR data. The first SFR declining for z ^ 1.5 is obtained from the original 
Durham parameters (see specifically Heyl et al. 1995). Their values a* = 1.5 and ah = 5.5 
minimize the star formation and maximize the feedback effects in the shallow, early wells; 
in particular, the value of corresponds to e D = v ~ 3 , that is, to galactic winds stronger 
in shallower wells. The corresponding L x — T relation at z — 0, 1 is shown in the middle 
panel of fig. 1 (left), and the prediction for the soft X-ray counts by the solid line in the 
bottom panel of fig. 1. 

The other extreme evolution is the flat SFR, also shown in fig. 1 (top); this is 
derived from the second set of parameter values listed in the caption. The set still leads 
to comparable optical luminosity functions but is more in tune with recent feedback 
reappraisals toward the low side (see Thornton et al. 1998, Ferrara & Tolstoy 1999; Martin 
1999); it includes the "neutral" values a* = and ah = 2. These correspond to the SF 
time r* oc v~ a * and the kinetic fraction e Q oc v 2 ~ ah being even for shallow and deep wells; in 
particular, the values of Vh, Q-h correspond to e G = const~ 10 _1 in the shallow wells, that is, 
to weaker winds. In the middle panel of fig. 1 (right) we show the corresponding Lx — T 
relation, and by the dashed line in the bottom panel of fig. 1 the corresponding counts. 

We stress that the deep counts in fig. 1 (bottom) draw increasingly apart for fluxes 
F x < 1CT 14 erg cm -2 s -1 ; in fact, at fluxes 10 times fainter the counts for weak exceed 
those for strong winds by a factor at least 3. Following up the discussion ending §3.2, 
the result can be reckoned on the basis of L x oc mlnl/nl (see eqs. 3 and 4). Here m h 
oc Z" 1 ; in addition, n\/n\ is enhanced at low T by about f~ x for weak compared to strong 
feedback; so the luminosities are enhanced as L x oc /~ 3 at low T. The numerical values of 
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/* differ by —20% at L x ~ 10 43 erg/s and In sum, with weak feedback we expect the 

deep counts, which rise here close to iV(> F x ) oc F^ 3 ^ 2 , to exceed the other case by about 
(0.8)~ 9 / 2 ~ 3 at Fx ~ 10~ 15 erg cm~ 2 s -1 , as confirmed by fig. 1 (bottom). 

5. Discussion and conclusions 

Having in mind surveys from Chandra and XMM we have limited our computations to 
photon energies exceeeding E min = 0.25 keV, and the resolvable sources to L x > 10 43 erg/s 
with Fx > 10^ 15 erg/cm 2 s. So we obtain a lower bound to the divergence of the counts; 
this is amplified by some 20% if luminosities down to Lx ~ 3 10 42 erg/s are included. In 
the critical universe, a given Lx would correspond to a brighter flux, to yield even larger 
excess counts at given flux. 

Considerable uncertainties currently plague the evolution of the SFR to z > 1, as 
derived from the optical (dust-obscured) and from the IR data (time-consuming to obtain, 
and with redshifts difficult to determine); these are discussed in detail, e.g., by Madau, 
Pozzetti & Dickinson (1998), Pettini et al. (1998), Hughes et al. (1998), Barger et al. 
(1998), Ellis (1998). So complementary information concerning the X-ray counterparts 
from the forthcoming surveys planned with Chandra and XMM will be welcome and timely 
to address the issue. Eventually, with the O-IR data consolidating in their own right, the 
soft X-rays will elicit the SN feedback in action. In fact, this contributes directly to the 
hot component, while it only indirectly affects the stars which form from the cooling stuff; 
therefore the soft X-rays are best suited to probe this key quantity. 

In fig. 2 we show not only the evolution of the stellar and the hot components from our 
computations, but also the abundances at various z of the third phase constituted by the 
lukewarm baryons inside and outside the galaxies. These turn out to be in fair agreement 
with the densities (actually lower bounds) evaluated from the observations of the Lya 
clouds seen in absorbtion, whether damped or not. 
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In conclusion, a weaker feedback AmJAra, allows more baryons (hot and cold) to be 
retained within shallow early wells, as illustrated in fig. 2. As shown in fig. 1, such baryons 
yield not only a larger star formation but also stronger X-ray emissions. These involve 
mainly outputs Lx ~ 10 43 erg/s or less at z « 1 or larger, associated with halo masses 
M ~ 5 10 13 M Q or smaller and with temperatures T ~ 1 keV or cooler. The deep counts in 
soft X-rays of souch sources and their Lx — T relation can probe directly the feedback, the 
main unknown that hampers the understanding the cosmic star formation history. 

We thank E. Giallongo, P. Rosati, and P. Tozzi for stimulating and informative 
discussions, and a referee for helpful comments. Work supported by grants from ASI and 
MURST. 
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FIGURE CAPTIONS 



Fig. 1.— Top: The z-dependent SFR from our model of baryon processing during the 
hierarchical clustering in the flat universe with Q Q = 0.3, H = 70 km/s Mpc. The solid line 
represents the shape peaked at z ~ 1.5 resulting from the Durham set of parameters: r ° = 2 
Gyr, a* = —1.5, Vh = 140 km/s, ah = 5.5, implying strong galactic winds. The dashed 
line represents the SFR flat toward high z resulting from our "neutral" set: r° = 2 Gyr, 
a* = 0, v h — 140 km/s, a h = 1.5, implying weak winds. Middle: The L x — T correlations at 
z = (solid line) and at z = 1 (dotted line) corresponding to the peaked (left frame) and to 
the flat SFR (right frame) illustrated on top. Group data from Ponman et al. (1996, solid 
squares); cluster data from Markevitch (1998, open circles). Bottom: The source counts 
iV(> Fx) corresponding to the peaked SFR (solid line) and to the flat SFR (dashed line), 
in the energy range 0.25-2 keV. The hatched band represents the cluster number counts 
observed by Rosati et al. (1998). 

Fig. 2.— The contribution to the baryon density Qf, from the lukewarm baryons with 
T < 5 10 4 K, from the hot virialized ICM and from the stars, see labels in the figure. The 
solid lines show the model predictions corresponding to strong winds, while the dashed lines 
correspond to weak winds. The density of lukewarm baryons is compared with the data 
(from Giallongo et al. 1997) concerning the density of photoionized intergalactic gas, which 
constitute lower limits especially for z < 3. 
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